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ABSTRACT 

A free-floating wave meter was used to obtain a time 
series of wave-induced pressure variations in deep water. 

The meter T s response was roughly checked against a piston- 
type pressure testing apparatus. 

Techniques in the operation of the meter were developed 
involving a slight modification in the suspension of the 
sensor . 

Power spectra of the pressure variations were obtained 
from the digitized wave record via the amplitude spectra 
given by the fast Fourier transform. 

Characteristics of the power spectra of the metered 
waves were used to evaluate the meter by comparing them with 
the waves that might be anticipated (on the basis of 
theory) from the existing surface wind field and with visual 
observations of the metered waves. 
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I. OBJECTIVES 



The U.S, Naval Postgraduate School acquired a free 
floating wave meter designed to measure deep water wind 
waves in the sea over the band of frequencies associated 
with swell and locally wind generated waves. It is desir- 
able to make operational use of an instrument of this kind. 
Most wave metering requires a physically stable reference 
platform from which to record deep water wave activity, 
or it requires that the reference platform's motion be 
known. These tedious and expensive requirements are avoided 
with the free floating meter. In order to facilitate the 
use of this instrument in future research, the objectives 
of this paper are: 

1. To estimate the response of the instrument by a 
direct calibration procedure; 

2. To providB a detailed account of techniques and 
equipment used in making wave measurements at sea with 
the instrument ; 

3. To analyze the wave records and to compare the 
results with independent estimates of the wave parameters 
made for the same wave system. 
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II . THEORY 



The free floating wave meter has a pressure sensor hung 
below a float which moves vertically with the free water 
surface . 

The wave pressure for wave number, k, of the fixed depth 
z^, is given by: 

-kz 1 

P w (k; z-^t) = p o (k;t)e 

where p Q is the static (or slowly time varying) component of 
the pressure due to the rise and fall of sea level. 

For the bobbing meter suspended at the nominal (mean) 
depth, , below the undisturbed sea surface, there are two 
components to the pressure variation: the dominant contri- 

bution is that due to the rise and fall of the sensor as the 
surface float is displaced, this being identically -p Q (k,t), 
and the other is the change in environmental pressure p . 
Thus, the total pressure change experienced by the meter is 
approximately: 

-kz-. 

P z (k;z 1 ,t) = -p Q (k;t) + p Q (k;t)e 
— k z 

p z( k;z i’t) = - P 0 U -e 1 ) 

Therefore, assuming that the sensor hangs vertically from 
the float, the instrument will detect pressure fluctuations 
corresponding to the vertical motion of the float. 

If one wants to measure the sea surface displacement as 
indicated by the bobbing motion, the meter must be suspended 
at a sufficient depth z^ , such that: 
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P z (k;z 15 t) = -p Q (k;t) 

That is, p z accurately describes the surface wave action if 

-kz -k z 

e 1 «1 for wave numbers k>_k Q . Letting e ° = .1 and 

z = 300 feet (the depth of suspension for the sensor in the 

measurement made in this study), pressure measurements from 

waves with wave numbers greater than .0077 feet~\ or 

frequencies greater than .079 seconds"'*' will not be in error 

by more than 10$. Gaul and Brown [1] give more detail in 

the theory of a free floating wave buoy. 
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III. DESCRIPTION OF THE INSTRUMENT SYSTEM 



The instrument system has a transducer-buoy package, 
a signal conditioning deck unit, and recording units 
(Figure 1 ) . 

The transducer-buoy package (Figure 2) consists of a 
sensitive pressure transducer, small diameter (.060 in.) 
sea cable, and buoy hull, which contains a modulator, 
radio transmitter, batteries, and antenna. The transducer 
detects pressure changes by means of a vibrating wire whose 
frequency varies directly with the applied hydrostatic 
pressure. The sensitivity and frequency response are 
discussed in Section IV. The transducer and buoy electronics 
convert frequency changes into an FM signal with a center 
frequency of approximately 27 MHz. The components of this 
package are designed to minimize water and wind drag. The 
transducer has a streamlined torpedo shape of length 30 cm. 
and diameter 3 cm. The buoy is circular with a 60 cm. 
diameter in horizontal section, and has an elliptically 
shaped vertical cross-section of depth 1*1 cm. 

The signal conditioning deck unit consists of a modified 
citizen’s band radio receiver, a discriminator, a d.c. 
amplifier, and a band pass amplifier. The deck unit condi- 
tions the signal from the buoy so that a ten foot peak-to- 
peak wave results in a five volt peak-to-peak signal. The 
output signal does not represent absolute pressures; it 
reflects changes from the static pressure. 
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Figure 1. Instrument System 




Figure 2. Transducer-buoy Package 



The recording unit provided with the instrument is a 
Leeds & Northrup Speedomax H chart recorder. In addition 
to the chart recorder, the output signal from the deck 
unit was sent through a d.c. amplifier and recorded on a 
Sangamo 3500 1^-track magnetic tape recorder. A signal 
generator was also used to put a 60 Hertz reference signal 
on a track other than the data track during the period of 
recording in order to identify the data to be processed. 

A double-beamed oscilloscope was linked to the playback 
outlets of the tape recorder in order to verify that the 
data and reference signal were being recorded. 

More detailed information about the instrument can be 
found in Gaul and Brown [2] . 
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IV. CHARACTERISTICS OF THE SYSTEM 



This section discusses four characteristics of the 
instrument: frequency measurement limitations due to size 

of float, calibration dial characteristics, frequency 
response of the transducer-buoy electronics, and the 
response factor of the instrument. 

A. FREQUENCY MEASUREMENT LIMITATIONS DUE TO SIZE OF FLOAT 
The accuracy of pressure measurements, in part, depends 
upon the ability of the buoy to follow the surface of a wave. 
The buoy does not respond to waves with a length less than 
twice the diameter of the float (2,33 feet). Therefore, let 
the minimum observable w ave length, L , equal 4.66 feet, 
and utilize the deep water equation 



where T is the period of the wave in seconds and L q is in 
feet. From this equation, T is 0.954 seconds or frequency 



not respond fully to passing waves with frequencies greater 



B. CALIBRATION DIAL CHARACTERISTICS 

The deck unit has a calibration dial used to maintain 
a 5 mvolt peak-to-peak output signal, which is equivalent 
to full scale deflection on the chart recorder for the same 
maximum wave height when the transducer is used at various 




is 1.05 sec ^ . Therefore, theoretically, the buoy does 



than 1.05 sec 1 (i.e. ~ 1 sec. waves). 
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depths. For example, a dial setting of 500 when the trans- 
ducer Is at a depth of 300 feet gives a 5 mvolt peak-to-peak 
signal when a ten foot peak-to-peak wave passes. Figure 3 
shows the dial setting for various depths to get full scale 
deflection for 10 foot waves. Figure 4 and 5 show dial 
settings vs. peak-to-peak voltages at constant depths for 
10 foot waves . 

C. FREQUENCY RESPONSE OF THE TRANSDUCER-BUOY ELECTRONICS 

Table 1 on frequency response of the transducer-buoy 

electronics was obtained from the manufacturer of the system 
Bissett-Berman Corporation, Hytech Marine Products. The 
table indicates that the pressure transducer is more respon- 
sive to pressure changes as the absolute pressure increases. 
By applying increments of pressure at various absolute 
pressures with a piston-type pressure gauge tester, it was 
concluded that the instrument does not record the magnitude 
of the pressure changes with sufficient sensitivity unless 
the absolute pressure is greater than 40 to 6 0 psi. There- 
fore, the instrument does not function properly in water 
with a depth less than 100 feet. This characteristic made 
it impossible to calibrate the meter by direct comparison 
with meters placed at near shore locations. 

D. RESPONSE FACTOR OF THE INSTRUMENT 

An attempt to compute the response factor of this instru 
ment was made by obtaining power spectra of a known input 
and its response. The piston-type pressure gauge tester was 
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Figure 4. Output Voltage vs. Dial Setting at 300 feet for a 10 foot 
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Figure 5. Output Voltage vs. Dial Setting at 150 feet for a 10 foot Wave 
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FREQUENCY RESPONSE OF TRANSDUCER-BUOY ELECTRONICS 



used to apply a 2 psi Impulse over 2 seconds. After the 
impulse, the pressure was held steady for approximately 20 
seconds . 

Three separate trials were made for which the input, as 
constructed from visual readings of the pressure gauge and 
stopwatch dials, is essentially as shown in Figure 6. The 
power spectrum analysis of this input and the analyses from 
the three outputs are compared in Figure 7. Variations in 
the spectra of these three responses were due to the dif- 
ficulties involved in creating three identical pressure 
impulses with the pressure gauge tester. 



Pressure 




Figure 6. Input from Pressure Tester. 
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The response factor was computed by the following 



formula : 



AREA. 
in 

AREA 



out 

R = Response Factor 

Area. = Area under input spectrum in stated frequency 
1 band 

Area Qut = Area under response spectrum in same band 

The results of this analysis were: 

1. For frequencies greater than .15 sec" 1 , the response 
factor was equal to 1.00. 

2. For frequencies between .10 sec" 1 and .15 sec" 1 , the 
response factor was equal to 1.21. 

3. For frequencies between .05 sec" 1 and .10 sec" 1 , the 
response factor was 1.63. 

Therefore, these results suggest that the instrument over- 
responds to low frequency (.05-. 10 sec" 1 ) waves. Frequencies 
less than .05 sec” 1 were not considered because the instru- 
ment wasn’t designed to measure frequencies that low. There 
is uncertainty in this analysis due to the uncertainty in 
the slope of the input pressure curve. 
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V. FIELD OPERATIONS 



This section deals with the problems and techniques 
involved in handling the instrument system at sea. 

The main difficulties at sea occurred in handling the 
sea cable. The sea cable has a tendency to coil and kink 
whenever it is not under tension. When a severe kink or 
puncture in the insulation occurs, the sea cable does not 
conduct the signal properly. Furthermore, the sea cable 
becomes difficult to handle when wet because of its small 
diameter. In addition to these problems, the 90 lb. test 
sea cable may snap with sudden jerks. 

These problems were solved by taping a small diameter 
steel cable foot by foot to the sea cable. The rigidity 
of the steel cable prevents the sea cable from coiling and 
the tape enables the operator to get a good grip on the 
line. Furthermore, the steel cable increases the safety 
margin against losing the transducer. 

The field operations described here were made from the 
U.S. Naval Postgraduate School's 6 3 foot oceanographic 
research boat. This boat is well-suited to operations 
involving this instrument because the launching and recover- 
ing of the buoy could be done by hand. A larger vessel 
would need a water level platform or smaller boat in order 
to handle the buoy easily. It would be difficult to use 
a hoist because of the delicate nature of the instrument 
and the large antenna mounted on top of the buoy. 
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In launching the floating package, the first step is 
to secure the buoy end of the safety line. With the boat 
engines stopped, the transducer and sea cable are then 
lowered in a hand-over-hand method, care being taken not to 
allow the sea cable to scrape against anything. This is 
best accomplished with the transducer lowered on the wind- 
ward side of the boat. Then the sea cable and safety line 
are attached to the buoy, which is lowered by hand into 
the water. To obtain a vertical alignment with the trans- 
ducer, the buoy is pulled away from the boat. 

An audible tone from the receiver indicates that the 
system is in operation. It is recommended that the boat 
maintain a position 30 to 60 yards downwind from the buoy in 
order to insure strong reception of the buoy signal and to 
keep the buoy in sight. After the boat is in position and 
the equipment is functioning properly, the reference signal 
generator is turned on to indicate which portion of the 
magnetic tape has data to be processed. 

The recovery of the buoy is similar to the launching 
operation. The buoy is easily brought aboard, by hand, with 
the base of the antenna serving as a handle. The sea cable 
and transducer are then hauled slowly aboard hand-over-hand. 

These techniques for handling the instrument were satis- 
factory for waves as high as 12 feet with white caps present 
For higher seas, different techniques in handling the instru 
ment may be required. 
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Further detail of the operation of the instrument can 



be found in the Bissett-Berman Manual [3]. 
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VI . RAW DATA 



Four wave records were made under varying sea and wind 
conditions. All four records were made in an area centered 

10 miles WNW of Monterey Harbor. All measurements were 
made in water with a depth in excess of 800 fathoms. Table 

11 describes the characteristics associated with each wave 
record . 

The Douglas Sea State Code taken from Fairb ridge [4] , 
describes codes 2, 3, and (used in Table II) as follows: 



State 2 small wavelets, crests of glassy appear- 

ance, not breaking; 

State 3 large wavelets, crests begin to break, 

scattered whitecaps; 

State k moderate waves taking longer form, many 

whitecaps, some spray. 
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DATA CHARACTERISTICS 



VII. DATA PROCESSING 



Processing the data consisted of digitizing the wave 
record and obtaining a power spectrum. 

Because of certain characteristics in the digitizing 
process, which will be explained later in this section, 
the data on the magnetic tape were played back at a speed 
of 6 0 ips (inches per second) (recording speed was 1 7/8 
ips) through a chart recorder moving at a speed of 
50 mm/sec. Furthermore, the data were played through a 
Khronhite filter, model 3321, with the cutoff frequency 
of 8 cps real time in order to get rid of noise introduced 
by electrical interference from the boat. 

The wave records (Figure 8) were digitized using a 
Calma Company Model 480 Digitizer which records Ay (of 
the cartesian co-ordinates of the pressure curves) for 
every .01 inch of- stylus travel in the x-direction. The 
digitized data were recorded on 5 inch mini-reels at 556 
bytes per inch. In order to make the data compatible to 
the IBM 360, the record on the mini-reel was converted to 
punched cards by program "Convert M (see Appendix B) on 
the CDC 6500 at the Fleet Numerical Weather Central, 
Monterey, California. Because of a limited buffer capacity 
program "Convert" is limited to 15000 sample points. There 
fore, 150 inches is the maximum wave record length that can 
be processed. This was the reason for condensing the wave 
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Figure 8. Filtered Wave Record (major division represents two feet) 



record by playing back the data at a speed increased by a 
factor of 32 times the recording speed. Under these con- 
ditions, one second of real data was represented by 1.56 mm 
on the chart paper. 

Due to the fixed sampling rate of the digitizer, the 
sampling interval. At, was .163 seconds. According to the 
Nyquist criterion, the maximum allowable At is: 



At 



max 



1 



2f 



max 



With f equal to 1.05 cps from Section IV, At = .^76 
max max 

seconds. Therefore, the actual sampling interval, .163 
seconds, successfully meets the Nyquist criterion. 

The power spectrum was obtained by utilizing IBM sub- 
routine RHARM. This subroutine computes one-dimensional 
Fourier coefficients by using the fast Fourier transform. 
The power spectrum was smoothed over by averaging every 20 
Fourier coefficients resulting in approximately *10 degrees 
of freedom, The entire procedure was performed by program 
’’Spectrum" in Appendix B, 
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VIII. COMPARISON OF RESULTS WITH THEORETICAL WAVES 



Figures 9-12 in Appendix A give the power spectra of 
the four records taken by the wave meter. These spectra 
were compared with visual observaitons and with waves 
expected from storm areas indicated by weather maps obtained 
from Fleet Numerical Weather Central. 

Significant height estimates were obtained by observing 
by eye (with two different observers) and recording the 
vertical movement of the buoy from trough to crest of a 
large number of successive waves. The average of the high- 
est one-third of these observations gives the significant 
height. From the power spectra, the significant height was 
calculated by 

H 1/3 = 2.83 VeT 

where 

E = 2 (Variance) = area under power spectrum 

Table III and Figure 13 give the results of these computa- 
tions. The results suggest that the meter may have responded 
accurately to both principle spectral bands. This is not 
in agreement with the analysis concerning the known input 
and its response from Section IV. However, Figure 13 is 
only suggestive because of the. limited number of records 
taken . 

The power spectra were compared to the reports from 
the Fleet Numerical Weather Central of the combined sea 
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Number 
of Waves 
Observed 


Visual 

H l/3 
( feet ) 


Variance 


Meter 

H l/3 
( feet ) 


May 11 


125 


5.10 


1.85 


3.85 


May 20 


75 


8.35 


9.15 


8 . 55 


J une 1 


288 


2.21 


1.45 


3.40 


J une 2 


261 


4.55 


2.03 


4.05 



TABLE III 

RESULTS FROM SIGNIFICANT HEIGHT ANALYSIS 




Figure 13. Eye Observations vs. Meter 
Observations of Significant Height 



3 ^ 



height, which gives contour lines representing significant 
sea height. 

All the combined sea height analyses for the appropriate 
days indicate a significant height of approximately three 
feet in the vicinity of the mouth of Monterey Bay (Figures 
1*1-17 in Appendix A) . This roughly corresponds to all the 
reports except the one taken on the 20th of May. However, 
the combined sea height analysis may not be too reliable 
because it is based on a small number of reports for the 
area along the coast of California. 

In addition, sea surface atmospheric pressure analyses 
were used to locate and investigate the generating areas 
of the dominant frequency bands for the swell recorded by 
the instrument. Taking the peak frequency from the power 
spectra, the group velocity of the dominant swell was com- 
puted using 

c = hl- 
g 2f 

Cg = deep water group velocity in knots 
f = peak frequency in sec"^. 

Then, the distance traveled by the dominant wave was computed 
by using the difference between the time of the sea surface 
pressure analysis and the time the wave record was taken. A 
line was drawn on each sea surface pressure analysis (Figures 
18-21 in Appendix A) to represent the starting location for 
the swell observed at Monterey Bay. These lines appear to 
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lie in each case within the possible generating area of the 
dominant frequency swell. 



Computations were also made to determine whether, accord- 
ing to a commonly used wave generation model, winds in these 
possible generating areas produced the recorded dominant 
frequencies. The geostrophic wind speed was calculated from 
the pressure field using 



V = geostrophic wind speed in knots 
S 

ft = angular velocity of earth 
<|> = latitude 

AP = pressure difference in millibars 

AN = latitude change ( 6 0 miles ^ 1 degree of 



latitude ) 

The surface wind speed, V g , is then approximated closely 
enough for present purposes by; 

V = ,8V 
s g 

With values of surface wind speed and estimated fetch, the 
maximum periods which these generating areas were capable 
of producing was estimated by using co-cumulative spectra 
for wind speeds as a function of fetch, as given by Pierson, 
Neumann, and James [51. In all four cases, the generating 
areas apparently were capable of producing the dominant 
frequency swell recorded by the meter. Table IV gives the 
results of the computations made. 




2ftpsin 



1 




AP = 9.33 AP 
AN sin<J> AN 



where 
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TABLE IV 



RESULTS FROM GENERATING AREA ANALYSIS 





Peak 
Freq . 

( sec~^) 


Period 

(sec) 


Travel 
Time 
(hrs . ) 


Miles 

Traveled 


V 

s 

(knots ) 


Fetch 
(miles ) 


May 11 


.08 


12.5 


4l 


900 


14.7 


200 


May 20 


.11 


9.1 


35 


550 


8.2 


200 


June 1 


.08 


12.5 


65 

fr 


19*10 


9.4 


200 


June 2 


.08 


12.5 


65 


19^0 


13.0 


200 



Furthermore, the surface pressure analysis seemed to 
indicate that secondary peaks in the power spectra are due 
to local wind conditions. For example, it is noted that a 
large secondary peak not present on 1 June occurs at 
approximately .25 sec ^ on 2 June. From Section VI, It was 
recorded that there was a large increase in the wind over 
that 24 hour period. This suggests that the secondary peak 
was due to the local winds. 
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IX. CONCLUSIONS 



The portable free floating wave buoy provides a rela- 
tively inexpensive method for gathering wave records. The 
evaluation of the records suggests that the main spectral 
bands, the swell and the local wind waves of the wave 
system, are properly resolved by the meter. However, 
further dynamic and static calibration is suggested in 
order to resolve the discrepancy in amplitude measurements 
between the input-response analysis and the significant 
height analysis. Comparison of the buoy's wave records with 
the simultaneously made records of well-calibrated wave 
meters on Flip would be a desirable future project. Also, 
further comparisons with weather analyses would be profit- 
able in evaluation of the meter. 

The portable nature and ease of operation of the free 
floating wave buoy provide the opportunity for deep water 
wave forecast verifications to be made by organizations 
which are financially limited. The meter permits records 
of deep water waves to be taken more quickly, easily, and 
cheaply than with most other existing deep water wave 
meters . 
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Figure 9. Spectral Wave Analysis ... 11 May 70 

90$ Confidence: = 1.5S, = . 71S 
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•05 .10 .15 .20 .25 .30 

frequency (sec”l) 

Figure 12. Spectral Wave Analysis ... 2 June 70 

90$ Confidence: = 1.5S, S g ^ = . 71S 




Figure • Combined Sea Height Analysis ... 
00Z, 12 May 70 

*3 




Figure 15. 



Combined Sea Height Analysis . . . 
00Z, 21 May 70 







Figure 16 . 



Combined Sea Height Analysis . . . 
00Z, 2 June 70 
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Figure 17. Combined Sea Height Analysis . . . 
00Z y 3 June 70 
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Figure 18. Sea Surface Pressure Analysis 
OOZ, 10 May 70 
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Figure 19. Sea Surface Pressure Analysis ... 
06Z, 19 May 70 
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Figure 20. Sea Surface Pressure Analysis ... 
00Z, 29 May 70 
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Figure 21. Sea Surface Pressure Analysis ... 
00Z 3 30 May 70 
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PROGRAM CONVERT 
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FLAG INSERTED IN "’’HE RECORD BY THE PERSON DIGITIZING 



o 






i 

i'. 



U-i 

K 

<s 

ZD 

or 

ii i 

a 

2 



ILi 

o 

< 

O- 

l/) 

> 



<r O 

co o 

►-H 0 rH — ■ 

mmQ 

— 

CD s: rC 

r- t'0 10 

a cr 

' O Oh oOh C O 

r— IfM- ■(— r-'Or— 'i — I HHr^rH 

_J l-CD o+ O + O + 

OOh< oOCooO oCOC occ 
COh25.^h I ^HOyi-C^ 
5: *~cl — || ii ii n H 

XOarCLL CXccoxcrcx® 

WCaLLDO^XOO-XOD Li 



o 


CM 










r-1 


rH 








X 


rH 


rH 














c 






if 


CJ 


o 


o 






if 


r-H |— 


1— 


»£) 




»>— . — 


if 


rH 






o 


IZZfM 


■W 


rH O 


O 


O 


r~* 


»- ► r O 


if 


o 


cd 


K 


r-4 


rH *— 1 «-<("■- 




o 








II II II U- 


X 


1- — 


— 


O 


c 


~ r-H r-H vl 


X 


xo 


XO 


CD 


X 


*- •* *t— 1 


fC' 


occ: o 


a: o 






— — • — r- 


r 


0 + 0 


+ © 


— 


o 


HtHMX*~ 


X 


X o 


X o 


O 


CD 




c 


— UJ 


K U 


o 




»Zh © 


CM 


cozz X 


XX X 


o 


— 


— — — ► LO 


*• 


iiO o h 


X o 1-0 


in 


rH 


- r r. LL 


t-HlX. 


o O X >4 Zvt DX«4'Z o 


»-! 


OpH 


ir<v ir\x<i- 


rHXMf 



OUhCOOL 1-0 IDO 

OJf 1 “’LU |l Zr,Or( II Xr~C II 

cXZ O XX ox 

a if n.h OC II ChC C ii l-r-H omLIC? 
XCM^ZUl—l-ZL'h'-Z + Sl-XI- II 

II— X— X X— XZ>X— X — x 

romU-C iLO — OC Li 0 — 0 il a O Li- O ' <3 

^yHUHODOLHOuSMOHCS 



ciUJOO rH rHrHr-n LI CHIfl 
|-~ XU- CMX CVCMCM— — O vL— CM 
OZ o I HI -2 \~r-i 

~ h-XK<<r— h-< 
2UZZZHCLZZG 
*— r Mcrazcnaro 
O’'Z' 0 'CC'Ch a: C^Z 

Q- Q CL U IX O LA/ CL U_ tO li : 



imn vo x 

oo o o 

rH CM H H 



00 0 

o o 



o 



LT.fOvO CrHCM 
r-» r-4 i-H QOC 

CM CM rH sOOvO 



o 



52 



Z 

* 

* 



s 

3 

K 

Eh 

O 

w 

Ph 

CO 

S 

< 

K 

O 

o 

PS 

cu 



a 

o 



X 

3: 



IS 




w 




z 










z 






a 














CO 


LL 




U- 




0 


*— 








<1 


O 




•- 




z 










X 






— 




< 


— - 








CO 


CO 




cr 




«- 


ll-' 










*— • 




o> 














u. 


CO 




0 




3 












X 




>fr 




ui 


LL 








— 


_J 




w 




c 


z 








H 


< 




CO 




•w- 


*-< 








3 


z 0 




► 






—1 








0 


c U-v 




— 




-J 




CD 






CO 


CO 




co 




<r 


z 


U' 






a? 


_J X 




<0 




>■ 


ZD 


X 






3 


<r (M 




0 




CO 


L- 


CD 






co 


CO * 




r 




u 


< 










K- <tf 








K 


O 


< 






O 


O X 




co 




z 




h- 






Z 


UJ t— 




•- 




1 — 1 


O 


<3 






<r 


a. 0 










K 


O 








</> Z 




CC 




CD 










CO 


ZIL 








z 


CO 


<1 




IT' — — 


*— 1 


JDJ 




0 




*— < 


X 






*-< cr ro 


to 


<T cr II 




<VJ 




3 


s 


LO 




► 00 


X 


Zh-O 








a 


w 


< 




► CO ff 


3 


00—0 




> 


CVJ 


z 


to 






II r-l r-< 


<3 


mUJw (\! 




z 


r-i 


<r 


LLI 


Cj 




U LL 


z 


COCO CG(V 




1 — 1 


•- 


co 


O 


UJ 




CO CO Of' 




zco 11 




► 


i — 1 


». 


ID 


(/) 




1 — »- »■ 




LL) 00 




— 


II 


— 


H- 


3 




z— - — ► 




ZXZCD 


— 


0 


M 


z 


1 — 1 — 


w 




—ZD— O'- 


<3 


— 1 — <L 0 


M 


CP 


►- 


w 


_IZ 




<T 


er -0 ►H 


CO 


CO z 


1 


r—i 


— 




Cl *• 


O 


H 


CL — 3UJ— 3U-. 


H 


O to 0 


O 


CP 


*- 


CO 




< 


< 


II (S. II LLI 


O 


u x a m 




w 


— 


h- 


•ct II 


U! 


C 


-J_J — _JLL 


UJ 


ZiOCDCD 


X 


X 


LiJ 


zz 


— ♦— > 


CO — 




<r *-uj ►w- 


a 


O II ll; 


1 


*- 


3 


►— < *- 


ITU *- 


(\J 


_J 


1 z — a 


co 


CD— cO 


< 


— — — 


■h 


C 3 


•-<> — 


1 0 


<r 


O-J — -JZ 




V)Zm ► 


— 


(VI vO O — 


0 . Ui 


r<t-i 


XZiOv 


z 




CO 



lL.<^a 'JTr-O OK 
h* CL CD * ^ — 'rJir — LL 

D 3 U l jvu.Chw O 

a. <i o <jjw^ 

Zcoiuh-* - u/ujku a. co 

o coco cocr — Cccli 

oco<oo <t-zzc<m 

SOW t— — — 3 

<t 1 1 — Zt_><f OCtOt/Mr H Z 

CL Cl ZD — — *■ * 

cd au(/;_j_i_JiiucFi: 

OCh-U>Q<<j 5 '<£v< 
Ol2ZX^5'U-U H<H-U C- LI 
a. <CUJ o— COCOOOCOU-CO 
h- 



LOK 
~ < 

Cv 
<co 
u c 

CO LL 



CD 

UJ — — 
_J — CV,U 
<1 cX 
Z — to > 
hOU-K- 
CDC'^l/ 

— nT r-*< 

CO *'—_j 
O^K 
— <au 
CCIK 
<r<r cO 
LliUJ OCO 
a'au -3 



to *■ 

X 

CD *- 
IT\» 
sf 



O *-4 

Z+ (X 

»— o' 



II — K 

o^xao 

O’— 0 11 ZDO 

-hOO-Z 
U. <r II Cmt-K 

h? — in + h z 

— Z — 
aC-C — Cco 
SILXCXUO- 



-D 

3 



II 



S - X- U- 

Z -W ~ 
XXX 
<—10 — 10 — U 

o * c -o ► 
o- c - c - 

>dox- a 



h- 

< 

o 



X~ sT —VL>— ZDUJ 

COO — — h- — K— b- ZK 
vcr- I U CLL <LL <1 a\ — 3 
ZO -K-DK-Zh^EZ + hQ. 

1 — « II ”3 II ►— <C0*— Q — CO” 3ZZ 

II -3 0 II Z cooroCCLC.' II CC' 
->-D 00 ^ 3 :a 5 :u 2 ii.D(ju 



CO 


O 


O 


O 


O 


0 


O 


O 


0 


0 


0 


O 


O 


O 


in 


0 


0 


O 


0 


0 


IL 




O 




vK 


in 


0 


X 


00 


CL 



53 





_J 




X 


z 








<1 


to- 


H 










> 


H 


o 


o 








a 


CO 


z 


a 








UJ 


c 


UJ 










K 




_J 


o 








Z 


to 


- 


m 








•— 1 


H 


— 


<y 










to 


X 


LL 








ID 


a 












z 


— 


r-i 


» 








a 


co 


to 


rvj 








2 


•> 


_ 


» 








C 


O 


00 










to 


— 


CL 


u. 








to 


o 


o 


O 








•- 


to 


to 










X 


o 


- 


to 








IT' 


► 


X 










»- 


o 


in 


<r 








to 


— 


H 


> 








to 


in 


— 


or 








o 


o 


X 


LL 








z 


0 


X 


H 








o 


o 


X 


Z 








LJ — 


— 


— 


►— « 








min 


in 


T 










IH — 


o 


Z 


or 








m r- 


o 


or 


ll; 








— to 


O 


H 


> 






>- 


in n 


-to 


O 


— O 






cc 


oto 


LL G 


LL) 


LL 






<1 


ir.H 


-JLU 


CL 


Z O 






tO 


hZ 


HO r 


I/O 


► z 






to 


LL*-h 


►hIL 




f-< *-• 






LU 


- O 


H 


_J 


II z 






o 


ii a 


►-<0 


<r 


*- z 






LU 




►z 


z 


z 






Z 


x<r 


-J<r 


o 


— 00 uo 








hh — 


LL* 


—H 


— o 






to 


co<r or 


CD tO 


to 


HvO>- 






< 


20 ar >- 


<UJ 


z 


to^r-flD 








LU- LU <r 




LU 


CCLU 






tom 


_J - LL Ol 


-e> 


Z 


-O Z 








X ^ OT 


oo; 


to-H 


— r~J — ) 






-JO' 


zoir - < 


— LU 


o 


►-r -a 








-or •- to 


oz 




~XH 






U-C 


JC- - >- 


-u 


ID 


or>i o^ 






H 


LLGtOZ > l_> 


CD 


z 


LL * UJ 






o 


C'LUCDZ _ Z 


-o 


o 


'-COO. 






a: O 


-O' lif ►"-» LLJ 


arm 


to 


— oto 






UJOr-J 


h- Ch ► X (MO 


ILZ 


V- 


SO 




-j 


M C 


-OC»X O- — ‘K- C7LU 


-a 


X 


HU-- 




2 CO 


— r~l 


— XLL LJ - Z <J> •}( h- li l Z 


OO— 


■vU — 


— LUX 




-Z 


ZCMOO 


GlTtOCLLU - X CM —41- CC • 


C-’ ll Ch xg cl I — 




HD 


OCT - 


JO -- l/'- ^ 4- * — — U.,-h 


intro t-vinn 




Zto X 


t— < hZO 


LUO- * x II O W U l-H +LL II 


H — Zcn-ixm „x 




II + X 


K 00 II o 


ChCu iua r-.*-. * X — HZZLL'h 


X r - 


a -r~<XH 




CD 


O c»— OUJ 


* — oi- < i -it- -K ^uj^-icr + >-~h 


U-mcor -r-*- 


-mo 


o 


H-iur 


LU LLi |l 3 


^sCwo<izor^ — ^ rr— * r LL'ccxr o 


xDa h 'C 


i-UvO^O 


c 


O-ZD 


^CO-ZCM^-hHJQLii 0+0 V |- ^LU»--}1 II II JOri 


1 II zo — 


H U 


i^HZ 


HO 


mcozto 


CLT ' 0 


1 u. <i ll z -r. tl _ ,| || w , xcm , 


— »-* HU 


<cmuj<ri/) 


4- ll 


ii *-« ii 


l/OZO +h-rH 


IZH5 -G-JCM—t* UJ + UJ + ^(_^^ ( VU.U.UH*“ 


►-H_)ZH ZXHZ 


cvxcr 


-.cci 


mZCC 


1 Ol — — 1 t II CV II -0 II MjZHr* II ZX_J 11 




•or n. 


i— o: Z u> n ii zcm zZ“? 


XU. OwO II 


ii or c it <r<aa ii o ii o ii wo-^wu — —<100 D'C.<aorc4! aco n 


•-JD ii era- 


1- HDU UZh^a UUZO-)X^>-JcQUCOcCZcE£I:OCXU.UUO 2U -n 


ZLL zxzz toHOtnocD 




r— 1 












c 


in o O 


o 


o 


o 




o c 


in 


h o o 


o 


o 


in 




O in 


O' 


O' C •-* 


CM 


m 


CO 




vj- 




fH r— ( 


«—4 


r— < 


f—l 




r— < *— 1 


o 


o o 


o 




o 







5 ^ 



a 

o 



o' 

K 

O 



C- 

c — ■ 



or.H 









l-H 


<rc 




w 










O - 


zo 




Cl 










U»_J 


C 




— 










O- UJ 


1— Cl 






— * 








t r>co 


OIL! 




►— < 


H 








< 


ZQ 




w 


— 


o 






X-J 


U i* 




cc 


to 


o 






1 - - 


Z - 




- 




to 






cc 


»-*x 




— « 


- 


•— • 






c - 


OlO 




»—< 




► 






2:0 


— 


0 


*— • 


X 


o 




M 


</) ► 


u - 


0 


or 


O 


o 




•fc 


cc. 


zo 


0 


LC — - 


— - 


If' 




H-f 


-h- - 


ou> 


>0 


— IT 


II 




X 


X 


oa 


• to 


x-x 


w 0 


Ll.' 


— 


< 


•f 


-JU 


fV 


<— 


0 


O 


o 


s: 




a - 


XM 




C\J 


<1 


o 


*o 


to 


0— 


lO*- 


"O'- 


— a 


CL 


4 


► 


0 


DC Cff -Jf 


-CU 


cc 


to 


H 


OJ 


O 




►X 


0 LL 


Cl - 






II 


+ h 


wr- 


•r-K 


II 


cox 


C 


— • 


1 


Ox 


3Sn 


-XO 


►—■X 


c-'O 


O* 


U. Uj 




CU u 


<< ►H C 


C- U- 


1 — «-H 




zz 


O 


O XX 


: cc a. vC 


U 0 


C o_ 


si' 


r*-» C"j 



O 

oHM 

C>++ -) I ZOMV II ZCG>-I- JOlOrHZ-h 

II — I — II f»—s£:? I — »—» LI <1 - II f"- II » — ' U ' <t 

coz- 2 f + s: — «-»— kk^: 

s: ii n -W'-z n it -zjj^axr-* *-*z~or 
X*--J li 2 lLr-iC?t— U- orO<J<3 D'Ci'-t w-c -Co: CT 
i/)ZZ-}")m(jC'Xxu uouSLiHdcauria 



e 

o 



o 

c 



o 

o 

0 - 




55 



SUBROUTINE RHARM 



Given 2N real numbers, X , X. ... X 0 ,. ,, the subroutine 

o 1 2N-1* 

computes Fourier coefficients, a Q , a 1 , b^, a 2 , b 2 ,... 
a N — X 5 b N-l 5 a n bbe e Q uat i° n 



Xj - | a 0 + ^a k cos(^) + b k sin(^) 



+ | a N ( “ 1)J 



where 



0 , 1 , 



2N-1 



Therefore, with N inputs from a time series, output 
contains : 
a Q 

~ 2~ 3 b 0 = a l 5 a 2 3 b 2 3 *** a N/2 3 b N = 0 



to get the power spectrum. Program Spectrum gives the spec- 
traX estimate, S, equal to: 



|+i |+i 

X S = X 
N=1 N k=2 



a S + b k 



where T = record length. 
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